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In view of the measured Higgs mass of 125 GeV, the perturbative renormalization group evolution 
of the Standard Model suggests that our Higgs vacuum might not be stable. We connect the usual 
perturbative approach and the functional renormalization group which allows for a straightforward 
inclusion of higher-dimensional operators in the presence of an ultraviolet cutoff. In the latter 
framework we study vacuum stability in the presence of higher-dimensional operators. We find that 
their presence can have a sizable influence on the maximum ultraviolet scale of the Standard Model 
and the existence of instabilities. Finally, we discuss how such operators can be generated in specific 
models and study the relation between the instability scale of the potential and the scale of new 
physics required to avoid instabilities. 
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I. INTRODUCTION 

The recent discovery of a light, to current experimental precision elementary Higgs boson has a major effect 
on the entire field of particle physics. A narrow Higgs boson with a mass of 125 GeV allows for an extrapolation of 
all fundamental interactions, except gravity, toward higher energy scales in terms of perturbative gauge theories. 
In particular, the renormalization group (RG) provides a link between the electroweak scale and more fundamental, 
higher energy scales. Renormalization group considerations for the Higgs sector typically provide two types of 
limits. An upper limit on the Higgs mass arises from the Higgs self-coupling becoming large (triviality bound) [3], 
while a lower limit arises from stability considerations of the Higgs potential nig. The observed low Higgs mass 
suggests that the triviality bound will not be a problem for energy scales below the Planck scale. When it comes 
to the stability of the Higgs potential the situation is less clear [BHIO]. 

As for any potential we refer to the Higgs potential as stable if it is bounded from below. We can attempt 
to link this condition to the sign of the quartic Higgs self-coupling evaluated at different energy scales [1], but 
the details of such a connection need to be tested carefully. Throughout this paper we will refer to this link as 
the perturbative or perturbatively renormalizable approach, because it typically relies on the assumption that all 
operators in the Higgs Lagrangian have at most mass dimension four. The great advantage of this approach is 
that the running couplings are known to high precision in the perturbative Standard Model, so we can compare 
its precise predictions to the measurements in the Higgs sector. 

Alternatively, we can ask the question which additional operators, for example higher powers in the Higgs field, 
can render the Higgs vacuum stable and how such effective operators can be linked to particle physics models at 
the corresponding energy scales. Whereas the canonical irrelevance of higher-dimensional operators for low-energy 
observables is natural from a Wilsonian viewpoint, the presence of such operators is expected at high energies, 
once we approach the scales of an underlying theory. Perturbatively, such operators can for example appear as 
threshold effects from integrating out massive fields. Alternatively, they can arise from other non-trivial effects of 
a UV-completion. As such higher-dimensional operators do not appear in the usual, perturbatively renormalizable 
Higgs Lagrangian, we refer to this generalized approach as non-perturbative. Its main advantage is that it makes 
fewer assumptions about the Lagrangian at energies significantly above our experimental reach. 

Based on the measured Higgs mass and top quark mass the perturbative analysis suggests that the electroweak 
Higgs vacuum might indeed not be stable. In the perturbative analysis the instability of the Higgs vacuum can 
be linked to the energy scale at which the Higgs quartic coupling turns negative. At the three-loop level this 
scale comes out around I0^° GeV, where the precise value depends sensitively on the Standard Model parameters. 
Several physical interpretations of this scale are possible: first, we can require that the Higgs potential be stable at 
all scales. New degrees of freedom would then have to appear below or at 10^° GeV, changing the renormalization 
group evolution, and rendering the potential stable ini0- We review one such approach based on a Higgs portal 
with a scalar dark matter candidate in the Appendix mm- Alternatively, we can require that the absolute 
value of the (negative) Higgs quartic coupling remains small. New physics only has to set in at the Planck scale, 
tunneling rates to the true vacuum have to be small, and our electroweak vacuum is meta-stable [5]. The third, less 
studied option [T!^UTH] includes higher-dimensional operators which stabilize the Higgs vacuum beyond 10^° GeV. 
Such operators will appear in theories with physical cutoff or matching scales. We explore this third option. 

The higher-dimensional operators can be linked to particle physics models, including heavy new particles m- 
Again, we choose additional scalars coupled through a Higgs portal. Unlike TeV-scale new scalars, heavy particles 
do not rely on their renormalization group running over a wide range of scales to affect the vacuum structure. 
Instead, they change the full set of running Higgs self-couplings at high scales to modify physics at intermediate 
scales. 

Our starting point is a field theory with higher-dimensional operators but including only Standard Model fields. 
This theory is perturbatively non-renormalizable. Accordingly, there exists a finite UV-cutoff A in the general spirit 
of the Standard Model as an effective theory. An appropriate tool for including all quantum fluctuations in the 
presence of higher-dimensional operators as well as a finite ultraviolet (UV)-cutoff is the functional renormalization 
group m- Our discussion of the stability of the Higgs potential proceeds in three steps: 

1. To quantitatively study the stabilizing effects of higher-dimensional operators we need a model which reflects 
all essential features of the Standard Model. We construct and test such a model in Secs. Ill A1 and III Bl 

2. In this toy model we describe the stability conditions, compute the possible Higgs mass range, and analyze 
its fixed-point structure (Secs. iH^i^ - 

3. Finally, in Sec. |HI| we study explicit models with additional heavy scalars and fermions and determine 
under which conditions they stabilize the Higgs vacuum, while yielding Higgs masses below the conventional 
stability bound. 
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In Appendix [X| we discuss the effect of light states, in Appendix [B] we summarize the computation of tunnelling 
rates, and in Appendices [C| and [D| we review the impact of higher-dimensional operators and give a detailed link 
between the perturbative and non-perturbative approach to the Higgs-top renormalization group. 


II. GAUGED HIGGS-TOP MODEL 

The aim of this paper is to investigate Higgs mass bounds and vacuum stability in the presence of higher¬ 
dimensional operators and a finite ultraviolet cutoff. In addition to the standard perturbative running, we use 
functional renormalization group (FRG) methods as a tool to compute the running of an extended set of operators. 
Let us first set up a toy model that allows us to study the essential features of the Standard Model in the context 
of vacuum stability. 

As a starting point, we briefly recapitulate the main features of the Standard Model at one-loop level. In 
the introductory sections we use H for the actual Higgs scalar, while (f denotes a general real scalar field which 
can play the role of the Higgs field H in our toy model. Once we arrive at our toy model which quantitatively 
reproduces the Standard Model, we will again use H for the corresponding scalar Higgs field. 


A. Standard Model running 


The perturbative approach starts from the usual Higgs potential, generalized to an effective potential by allowing 
for a dependence on the momentum scale k of all parameters. Including higher-dimensional operators it reads 




A 2 n(fc) 
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where we denote the scalar doublet as {H/\/2,0) = {{h + v)/'/2,0), with h parametrizing the excitation above 
the expectation value v and neglecting the Goldstone modes. In this form it is not at all clear whether a slightly 
negative A 4 will lead to an unstable potential. This depends on the higher-dimensional couplings Ae^s,... which, if 
sufficiently large, can obviously stabilize the Higgs potential for all k < Mp\. 

Standard power counting shows that long-range observables are dominated by perturbatively renormalizable 
operators. Hence, one might think that higher-dimensional operators can generically be ignored from the outset. 
With the so-defined Standard Model the corresponding Lagrangian consists of all dimension-4 operators. In that 
case the question of stability is usually linked to the sign of A 4 defined in Eq.Q. The beta function for any 
coupling g is defined as /3g = dg/dlogk. In these conventions the one-loop renormalization group equations for 
the Higgs self-coupling A 4 , the top Yukawa y, and the strong coupling gg in the Standard Model read 
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The top Yukawa coupling is linked to the running top mass by y = V2mtlv, while the Higgs mass is given by 
= 2 A 4 U^ plus contributions from higher-dimensional operators. The couplings or the related running masses 
mt and rrih can be translated into the corresponding pole or on-shell masses, for example to describe the kinematics 
of production and decay processes. The gauge couplings are gs for SU{3)c, g 2 for SU{2)l and gi for {7(l)y. The 
number of fermions contributing to the running of the strong coupling is Uf. In the Standard Model setup of 
Eq. (§ no explicit higher-dimensional operators occur. However, if they are generated by Standard Model particle 
loops, the effect of, e.g., an induced As coupling is included in the running of the renormalizable parameters at a 


higher loop order. We will return to a more precise discussion of this point in Sect. HE 


To allow for an easy comparison of our toy model to the Standard Model we can translate the beta function 
for the Higgs self-coupling into the running Higgs mass, namely 
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Let us take a brief look at the essential features that arise from the renormalization group running given by 
Eq.([^. Toward the UV, the negative contribution of the top Yukawa term drives the Higgs self-coupling A 4 to 
small values and eventually through zero. Ignoring higher-dimensional operators, the Higgs potential seems to 
become unstable at large energy scales, before the stabilizing effect of the weak gauge coupling sets in at very 
high scales. This is the usual (meta-)stability issue in the perturbative setting. We will give a more detailed 
interpretation in Sect. ra 

For the quantitative behavior of A 4 in Eq. it is important that the top Yukawa coupling y also decreases toward 
the ultraviolet under the influence of the strong coupling as- Indeed, for the situation we want to investigate, A 4 
and y both decrease toward the ultraviolet, and the loop contributions from those couplings to their running get 
smaller and smaller. At some point the contributions from the weak couplings and y to the running of A 4 become 
comparable in size. When their loop contributions almost cancel, the running of A 4 becomes flat, j3\^ « 0. Finally, 
in the deep ultraviolet the weak gauge couplings dominate and turn A 4 back to positive values. This behavior is 
shown in the upper left panel of Fig. 

Following this argument, the self-coupling A 4 in the Standard Model appears to first turn negative and finally 
become positive again. From this behavior one could conclude that we live in a so-called meta-stable vacuum if 
the Higgs potential, including only operators to dimension-4 is given by 

TT4: 

H(iL) « A 4 (iJ) — with H ^ , (4) 

at high field values. Here, it is assumed that the effective potential is well approximated by identifying the RG 
scale with the field amplitude, X 4 {H) = A 4 (fc = H). This potential features two minima, our electroweak minimum 
and a global minimum at very large field values H ^ Mpi. The latter occurs far outside the region where the 
renormalization group equations of Eq.([^ can be trusted. We therefore focus on stabilizing effects that set in 
below the Planck scale. 


B. Toy model 


In this section we set up a simple model that exhibits most of the essential features of the behavior of the 
running Standard Model couplings without inheriting the full gauge structure. First, we replace the Higgs field H 
as part of an SU{2) doublet by a general real scalar held (p featuring a discrete Z 2 chiral symmetry. This ensures 
that no Goldstone modes alter the renormalization group flow in the symmetry-broken regime, just as in the 
Standard Model with the full gauge structure. In a simple Yukawa system without gauge degrees of freedom we 
already observe a similar perturbative flow toward negative A 4 at high scales M- The running of the top Yukawa 
coupling to smaller values in the ultraviolet is included when we add an SU{3) gauge sector. Gorrespondingly, we 
investigate the Euclidean action defined at the UV-cutoff scale A 


^A = 
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with an SU{Nc) gauge field {N,, = 3), a real scalar </?, Dirac fermions ipj for Uf flavors, the covariant derivative 
including the strong coupling gs, and the effective scalar potential I 4 ff- To avoid confusion, we emphasize 
that the Higgs potential I4ff(A) is the bare potential from the Standard Model point of view. Embedding the 
Standard Model into a more fundamental theory, the Higgs potential becomes the effective potential of that 
theory, containing all effects of quantum fluctuations above A. We assume that Uy < nj of the fermion species 
are heavy and shall have a degenerate large Yukawa coupling y. The remaining Uf — Uy flavors have negligible 
Yukawa couplings. The effective scalar potential is now expanded in terms of dimensionless couplings A 2 n as 


V,s{k) 
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This potential should be compared to the Standard Model potential shown in Eq. Q . As part of the potential we 
study beta functions for the self-interaction A 4 and the (/ 5 ®-coupling Ae. Here we restrict ourselves to the first in 
an (infinite) series of possible higher-order couplings. 


Because our toy model does not reflect the weak gauge structure of the Standard Model, the gauge couplings 
gi and 52 do not appear. However, we know that in the Standard Model they have significant effects: first, the 
gauge couplings give a significant positive contribution to , balancing the negative top Yukawa terms for small 
values of A 4 ; second, they decrease the Yukawa coupling in the UV, thereby also increasing A 4 toward large scales. 
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Figure 1: Upper left: running of SM conplings, figure taken from Ref. [^. Upper right to lower right: running of our toy 
model couplings after including the rnnning strong coupling and electroweak coupling effects. Dashed lines indicate the 
regime where A 4 (fc) < 0, which in the perturbative approach dehnes the loss of vacuum stability. 


Since the variation of the weak coupling at large energy scales is modest, we account for its effects by including 
a finite contribution in the beta functions for A 4 and y, parametrized by a fiducial coupling gp and numerical 
constants c\,Cy. 
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The scalar mass is given by = 2 A 4 W^, and the running top mass by y = \f2mtlv. To approximate the Standard 
Model we choose n/ = 6 to account for the contribution of all flavors to the running of the strong gauge coupling. 
As long as we only keep the top quark contribution to the running of the Higgs quartic coupling, we set Uy = 1. 
The expressions in Eq.Q reproduce the standard one-loop beta functions for our model if higher-dimensional 
operators are ignored jl9[ . 

We have verified that the simple modelling of the electroweak gauge boson effects in Eq. Q already provides 
a reasonable estimate. In a more self-consistent treatment, electroweak contributions to all running couplings 
/r^, A 4 , Ag ... of the Higgs potential appear. This is discussed in more detail in Appendix 


The running of the strong coupling and the top Yukawa coupling in our model agree exactly with that of the 
Standard Model given in Eq.([^. Only the running of the scalar quartic vs Higgs quartic coupling is slightly 
different because of different numbers of degrees of freedom. To compare our toy model with the Standard Model 
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we also write the beta function for the scalar quartic coupling in terms of the scalar mass as in Eq.([^, 

3 

+ 2 TO(TO^ — 4m( + weak and Ag terms 
—- r(0.93m„)'^ + 2mirn? — Amf + weak and Ag terms] . 




( 8 ) 


The only difference is the self-interaction contribution to the running of the scalar or Higgs mass. This means 
that we can model the running of the top-Higgs sector if we slightly shift the scalar and top mass ratio at the 
per-cent level. 


In Figj^we illustrate our toy model with representative one-loop renormalization group flows in three steps: 

1 . in the upper left panel we show the full 3-loop Standard Model running from Ref. [3] as reference; 

2 . in the upper right panel we include a scalar field and one heavy (top) quark, Uy = 1, without any contribu¬ 
tions from the gauge couplings (pure Z 2 model); 

3. in the lower left panel we add the effect of as on the running top Yukawa coupling; the strong coupling itself 
runs with nj = 6 active flavors; 

4. finally, in the lower right panel we also include a constant fiducial coupling with c\ = 9/16, Cy = 97/30 and 
qf = 0.57 to model the effect of the weak coupling. 

Since this choice of parameters allows us to quantitatively reproduce the behavior of the Higgs sector in the 
Standard Model, we refer to the corresponding scalar field given the appropriate model parameters as the ‘Higgs 
scalar’, i.e. ip —>■ H. For the flow trajectories we choose the running top mass such that mt{mt) = 164 GeV [S], and 
correspondingly the Higgs mass at nih = 2\4^{mt)v^ = 125 GeV. This means that we ignore the small difference 
between the running Higgs mass mhijnt) and the measured Higgs pole mass. The running of the strong coupling 
satisfies as{mz) = 0.1184 at mz = 91 GeV [9j. 

When we include the strong coupling in the bottom left panel of Fig. the top Yukawa coupling decreases 
toward the ultraviolet. This allows for lower Higgs mass bounds from vacuum stability compared to the pure Z 2 
model [14]. With the additional approximate effects from the weak coupling we also reproduce the back-bending 
toward positive A 4 in the deep ultraviolet, as it occurs in the Standard Model, shown in the bottom right panel 
of Fig. The weak coupling term Cy in the flow of the top Yukawa coupling is responsible for the faster decrease 
of the Yukawa coupling, such that it flows almost parallel to the strong gauge coupling. Thereby also the flow of 
A 4 becomes flatter. The weak coupling contribution c\ to the running of the Higgs quartic coupling has a smaller 
impact at intermediate scales and induces the final back-bending of A 4 close to the Planck scale. 


C. The issue of vacuum stability 

In this section we clarify some aspects of the relation between the perturbative and non-perturbative approaches. 
First we address a potential confusion whether the top loop stabilizes or de-stabilizes the potential. Second, in 
the perturbative approach the RG scale is replaced by the field value. In the non-perturbative approach both, 
the RG scale and the field value appear explicitly. Finally, in the perturbative approach one usually employs a 
massless regulator based on dimensional regularization. The question arises what happens in the presence of a 
finite, physical UV-cutoff. 

Let us start by briefly recalling how the issue of vacuum stability arises and what role a finite cutoff plays. 
We refer to a potential as stable if it is bounded from below, otherwise it is unstable. A stable potential can 
exhibit several minima, where a local minimum might be meta-stable while the global minimum is stable. As 
argued above, the top loop is by far the dominant factor in the renormalization group evolution of the Higgs 
potential. Indeed, the whole issue of stability can already be understood qualitatively by considering the term 
/ 3 a 4 = - 3 ^/( 8712 ) -h • • • in Eqs.^ and 0 . Integrating the /3-function toward the ultraviolet, this term drives 
the quartic self-coupling to negative values at some high scale. This is the usual argument for the loss of stability 
in the Higgs potential if the Higgs potential is approximated by V « X 4 {H)H'^/A. 

This observation appears to be in conflict with non-perturbative lattice simulations |20] and with arguments 
that the interacting part of the top loop is non-negative and cannot induce any instability in the presence of a 
finite UV-cutoff scale A mi EH). To understand this in more detail, let us consider a one-loop calculation of the 
effective potential at the electroweak scale k « /cew ~ 0, where we neglect A:ew ^ A when appropriate. As 






7 


shown in Ref. |15j the fermion determinant representing the top contribution to the effective potential for any 
finite cutoff A can be written as 


AVtop = + positive terms, 


(9) 


where C 2 > 0. We start from a stable bare potential Vuv = f4ff(A) of the kind shown in Eq.Q with quadratic 
and quartic Higgs terms. The top loop then gives us an approximate weak-scale potential 


f4ff(fcEw) ~ V\]y + AVtop = 
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and will never lead to an unstable effective potential at the weak scale. A closer look at the contribution from the 
top loop below the cutoff scale A reveals the formally sub-leading terms 
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where we normalize the contribution such that AVtop(77 = 0) = 0. In the last line we adopt the general notation 
of Eq.®. The additional logarithmic contribution to the quartic potential term contributes to the Higgs mass 

as [Tilfe] 


A4(fcEw) = A4(A) -t C47/4log 


A 


1 


fcEW 8 


( 12 ) 


This is nothing but the renormalization of the Higgs self-coupling. For a sufficiently large cutoff scale A the above 
top contribution forces us to choose a negative high-scale value A 4 (A) < 0 to obtain the measured Higgs mass. 

While no instability in the effective potential is induced by the top loop, its contribution is such that we would 
have to choose an unstable potential at the ultraviolet scale A to reproduce the measured Higgs mass in the 
infrared. The usual perturbative renormalization group evolution simply determines the value we would have to 
choose in the ultraviolet in the class of pure iJ^-potentials and in that sense reflects the fact that we would have 
to choose an unstable UV-potential if we want to start at a very high scale A. 


Stabilizing the UV-potential requires higher-dimensional operators. In this situation the stability of the UV- 
potential cannot be determined from the running of A 4 alone, but requires the full UV-potential. As a simple 
example, we include a Ag term as in Eq. 0 / 


Vuv = Ueff(A) 


/ 72 (A) 
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As we have seen, A 4 (A) is essentially fixed by the measured value of the Higgs mass in the infrared. For sufficiently 
large A it is negative. Similarly is fixed by requiring v = 246 GeV in the infrared. For the simple momentum- 
space cutoff regularization, it is typically large and positive, y^{A) ~ C>(0.01)A2. As expected, A 6 (A) is the only 




Figure 2: Left: Sketch of the possible shape of Uuv with y?{A) > 0 but A 4 (A) < 0. Potentials are stabilized by A6(A) 
smaller (red) and larger (bine) than A4(A)A2/(3 /i 2(A)), respectively. Right: Sketch of different possible shapes of the 
effective potential « 0) that may occnr depending on the size of A6(A), all corresponding to a stable UV-potential 

with a single minimum at 77 « 0 as the blue curve in the left panel. Note that the resolution of this sketch is not high 
enough to display the offset of the electroweak vacuum from 77 = 0, the corresponding two minima appear as one at 77 = 0. 
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Figure 3: Starting from a stable potential in the UV different flows to the IR potential with the electroweak minimum 
are possible. In the left panel we show a situation where the potential remains stable during the whole flow from k = A 
to the point where the electroweak minimum forms. This corresponds to region I in Fig. For smaller initial values of 
Ae corresponding to region II in Fig. we observe a meta-stable behavior of the scale dependent potential already at 
intermediate values of k, which then disappears. The electroweak symmetry breaking minima are not formed until later. 
We call this behavior pseudo-stable. 


new and free parameter. Its value is essentially undetermined by measurements in the infrared, because it is 
suppressed by the large scale A. 

Choosing A 6 (A) positive ensures that the UV-potential grows at large field values, making it a viable bare 
potential for a quantum field theory. In principle, this can be achieved with an arbitrarily small value of A 6 (A). 
Since /i^(A) > 0, the potential then has a qualitatively similar form to the red curve in the left panel of Fig. ^ 
with a local minimum at H = 0 and a global minimum at a large field value H ^ 0. Choosing 

A 6 (A) > A2(A) (14) 

removes the second minimum at 7 ^ 0 for negative A 4 (A). The only minimum appears at Ff = 0, as shown by 
the blue curve in Fig. In this paper we are mainly interested in this type of potentials. 


Nevertheless, a stable UV potential with a single minimum at iJ = 0 according to Eq.(14) does not necessarily 
imply that the effective potential only features a single minimum. In the course of the RG flow, all types of behavior 
with a potential bounded from below shown in Fig. do occur, depending on the choice of UV parameters. 

Qualitatively, this behavior can be understood from the unique-vacuum condition in Eq.(14), which can be 
applied to any scale k. While Xe{k) decreases rapidly for decreasing fc, A 4 grows at the same time and ultimately 
becomes positive. If A 4 turns positive first, then the scale dependent potential is always stable during the RG 
evolution in our approximation, cf. left panel in Fig. On the other hand, if Eq.(14) is violated first, a second 
minimum appears and the minimum at iF = 0 can become meta-stable for some range of RG scales, cf. right 
panel in Fig. However, a reliable (meta-)stability analysis of the full potential will require to go beyond the 
simple polynomial expansion of the effective potential around one minimum. 

For completeness let us note what happens when we continue to run to lower energy scales. Independently 
of the high-scale behavior discussed above, drops below zero near the electroweak scale. For these low 

scales A 4 is always positive and the interplay with the negative results in a finite vacuum expectation 

value V. In the presence of higher-dimensional operators, we rather generically find RG flows that interconnect 
a stable UV-potential with a minimum at iF = 0 with stable effective potentials in the infrared and a global 
electroweak minimum at H = v. If the flow passes through a finite regime where our polynomial approximation 
looks meta-stable, we refer to the scenario as pseudo-stable. In that case a full stability analysis requires a detailed 
non-perturbative analysis of the effective potential. 


Our discussion shows that we have to distinguish between the stability of the UV-potential where no quantum 
fluctuations below A are taken into account and that of the effective potential with all fluctuations included. In 
principle the UV potential and the effective potential can have quite different shapes as sketched in the left and 
right panels of Fig. In the approximation V ~ A 4 (FF)FF"^/ 4 , this difference cannot be accounted for properly: on 
the one hand the running quartic coupling defines the UV-potential as U(A) = A 4 (A)FF^/ 4 . On the other hand, 
the identification fc —>■ FF is assumed to be a good approximation for the effective potential. This means that 
one and the same function X 4 {k) determines the stability of the UV-potential and the stability of the effective 
potential. The presence of higher-dimensional operators influences both aspects: at the UV-scale, they can modify 
the UV-potential in a rather general way. Successively, they can contribute to the RG flow of the renormalizable 
operators for some range of scales. The functional renormalization group takes these aspects into account by 
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describing a scale-dependent effective potential T4ff(fc; H), depending on both the scale k as well as the held value, 
such that 


V{k = A-,H) = Vuv and V{k = 0;H) = Ks- (15) 

In this manner, the (meta-)stability properties of the effective potential can be followed in a scale-dependent 
manner. 


Finally, let us note another subtlety in the comparison between the perturbative and non-perturbative ap¬ 
proaches. The perturbative approach usually relies on a mass-independent regularization scheme. To avoid large 
threshold effects one should therefore stop the running of the couplings at the appropriate mass scales. Most of 
the relevant masses are proportional to H. This suggests to approximate 


FeS = Vik = H,H)^ 


(16) 


as is usually done in the perturbative approach. In Eq.(I2) we have already identihed the running Higgs self¬ 
coupling at the scale /cew- Using the identification k = H we essentially reproduce the first term in the square 
brackets of the second line of Eq. Let us check this explicitly. 
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However, 14^ = /A only holds for iJ z; ~ |/r(fcEw)| and A^ H. This means that in the presence of 

a finite, physical UV-cutoff A the perturbative running of X 4 {H) provides a good approximation of the effective 
potential only for field values H A. Indeed, using the perturbative running for a cutoff-free mass-independent 
regularization scheme can make the effective potential appear to be unstable beyond H A even when it is 
not [m |20]. Essentially, this results from a breakdown of the expansion in powers of H/A when going from 
the second to the third line in Eq.([T7|). For the use of the perturbative approach we therefore always have the 
condition that the held values — more generally any scale under consideration — are smaller than any physical 
cutoff. 


D. Higgs mass bounds 

When it comes to stability considerations, the main novel feature of our renormalization group hows is the 
presence of a second model parameter, Ae, as a representative of general higher-dimensional operators. In the 
absence of Xq the stability condition for the bare potential, A 4 (A) = 0 , at the so-dehned ultraviolet cutoff scale 
A constitutes the lowest viable choice for the quartic coupling. Positive higher-dimensional couplings in the 
ultraviolet can allow for negative values of the Higgs quartic coupling and yet a stable potential. A negative 
quartic coupling in the ultraviolet also lowers the infrared quartic coupling value, leading to a stable potential 
with a Higgs mass below the conventional lower bound obtained when neglecting Ag. 

Let us now consider two examples for the inhuence of Ag > 0 as shown in Fig. In the hrst case, corresponding 
to the left panel of the hgure we assume a relatively large UV-value of Ag. At hrst, the effect of this term 
is sizable and therefore the RG how deviates signihcantly from the usual perturbative scenario. However the 
canonical dimension of Ag implies that it quickly decreases toward the infrared, while A 4 increases. At some 
scale the combined renormalization group how reaches the point Ag « 0 and A 4 « 0. This setup coincides with 
the ultraviolet starting point of the perturbative approach neglecting Ag and marks the boundary of stability. 
The infrared regime to the left of the dashed line in the left panel, where A 4 again reaches positive values, is 
basically identical to the usual perturbative scenario. This means that due to higher-dimensional couplings our 
scenario features a signihcantly increased cutoff scale, while reaching the same values for the Higgs mass, quartic 
self-coupling, and top mass in the infrared. The naive upper limit on the validity of the model near 10^*^ GeV 
(dashed line) is replaced by a higher scale ^ 10^^ GeV, at the expense of a nontrivial running of higher-dimensional 























10 



Figure 4: Running couplings including A 4 and Ae, corresponding to Eq. 0 with rrih = 125 GeV. The bare potentials 
satisfy the stability condition of Eq.( |14[ ) and thus remain bound from below at all lower scales. Left: A 4 (A) = —0.058 and 
A6(A) = 2.0 and A = 5 • 10^^ GeV. Right: A4(A) = —0.050 and A6(A) = 0.8 and A = 10^® GeV. 


operators. In between these two scales an effective field theory description within the model is still possible, while 
true new physics would presumably have to set in beyond 10^^ GeV, in order to control the strong running of the 
scalar couplings. 

In the second scenario, shown in the right panel of Fig. we choose a smaller value of As > 0 and a much larger 
cutoff scale A. In spite of this, the UV-potential at A is still stable with a global minimum aX H = 0. As we have 
chosen a significantly larger cutoff scale A, Ag becomes smaller than required by the unique-vacuum condition 
Eq.( 14) before A 4 reaches positive values. This implies an intermediate pseudo-stable region. The existence of this 
intermediate region indicated by dashed lines is a characteristic of the pseudo-stable scenario mentioned above. In 
our simple polynomial expansion, the scale-dependent effective potential develops a second minimum and H = 0 
seems to be only a meta-stable minimum for intermediate RG scales. Our analysis based on the functional RG 
as outlined in App. going beyond the simple perturbative arguments indicates that this behavior should be 
assessed critically. We tentatively consider it an artifact of the polynomial expansion of the effective pot ential 
around one minimum. The functional RG provides us with a /3-function for the whole potential, Eq.(C3). We 
evaluate it in a polynomial expansion about an expansion point H = 0 at high scales and H = v near the weak 
scale. After integrating the flow we check whether the solution for the scale-dependent effective potential satisfies 
this /3-function. Whereas the accuracy of the solution is indeed at the level of machine precision near the expansion 
point, the accuracy depletes away from this point. The convergence radius of this expansion is typically of the 
order of rrih M- In particular, near the apparent new minima the polynomial expansion should not be trusted. 
Reliable statements about the stability of the effective potential in pseudo-stable scenarios require a dedicated 
analysis of the RG flow of the effective potential. 

An improved functional RG analysis according to App. 0 also clarifies the following apparent problem: in 
the infrared, the Ag coupling runs negative seemingly indicating a global instability. This is prevented by even 
higher-operators ~ iJ®, which are taken into account in our improved functional RG analysis that we employ to 
determine the stability regions shown in Fig. [^discussed below. All our explicit solutions for the scale-dependent 
effective potential are bounded from below. 


Let us now examine the effects of Ag on the stability bounds for the Higgs mass. To derive Higgs mass bounds 
we first focus on a stable UV-potential with a global minimum at vanishing field values. In the absence of higher¬ 
dimensional couplings, the lower Higgs mass bound can be determined by requiring the quartic coupling to vanish 
at the cutoff scale A. Higher-dimensional couplings at the scale A can take values of 0{1) in units of the cutoff 
scale. Indeed, this is the generic situation that one would expect when examining an effective field theory such as 
the Standard Model close to its cutoff scale. In the left panel of Fig. we first confirm that the observed Higgs 
mass around 125 GeV corresponds to the ultraviolet boundary condition A 4 = 0 at A « 10^° GeV. Shifting this 
boundary condition closer to the Planck scale and hence allowing for a fully stable potential would require a Higgs 
mass close to 130 GeV with a fixed top mass [ 8 j. We also show several choices of ultraviolet boundary conditions. 
Those which predict smaller physical Higgs masses than the choice Ai(A) = 0 allow us to increase the ultraviolet 
cutoff. According to Eq.([l4|) a viable choice for a UV-stable potential including higher-dimensional operators is 
for example A 4 = —0.05 and Xq = 0.8. The corresponding Higgs mass stays below the conventional lower bound 
for all values of A and indeed gives a stable potential at the Planck scale, with a dip into a pseudo-stable regime 
at intermediate scales. 

As described in detail in Sec. |H B| our calculation of the Higgs mass relies on a numerically convincing, but 
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Figure 5: Infrared value of the Higgs mass (left) and the relative shift measuring the departure from the conventional 
lower bound (right), as a function of the UV-cutoff A. We show different UV boundary conditions as indicated in the plot. 


nevertheless a toy model. Special care is required when translating the computed shift of the Higgs mass at fixed 
cutoff to the Standard Model. We nevertheless conclude that shifts at the level of 1 — 5% seem viable. To see 
this we study the relative shift in the Higgs mass between the perturbative Standard Model and including the 
dimension -6 self-coupling in Fig. It is given by 


Amfc (A) 


(AdA)#0) _ (A.(A)=0) 

'"‘h _((bi_ 

™(^.(A)= 0 ) 


(18) 


The dependence on A is induced by the choice of ultraviolet cutoff scale at which we define the boundary conditions 
for A 4 or Ag. Negative values imply that the Higgs mass resulting from finite ultraviolet values of A 4 and Ag lies 
below the conventional lower bound without exhibiting instabilities in the potential in the UV or at the weak 
scale. We emphasize, however, that large negative go along with pseudo-stable scenarios that deserve more 
detailed investigations. 

For all choices of UV-couplings we observe that | | decreases with increasing A. This can be understood from 

the following argument [ 22 j : including higher-dimensional couplings can allow us to extend the renormalization 
group flow toward the ultraviolet by several orders of magnitude, i.e. the cutoff scale is increased beyond the naive 
estimate determined by A 4 (A) = 0. This corresponds to shifting the curve mh{K) in the left panel of Fig. [^to the 
left by some amount A A. The cutoff dependence of TOft,(A) flattens toward higher cutoff scales, thus such a shift 
is less effective for large cutoffs, and correspondingly A^^ decreases in this regime. 

While the achieved shifts in the allowed Higgs masses seem rather small, we emphasize that with present 
uncertainties in the experimental input parameters absolute stability of the Higgs potential without Ag terms is 
disfavored at ^ 3cr [^. Our shift of the Higgs mass limit should be viewed as a shift in the central value which 
signihcantly reduces the tension between the measured data and the possibility that the Higgs potential is stable 
up to the Planck scale. 


So far, we have restricted ourselves to absolutely stable bare potentials with a minimum at vanishing field values. 
Next, we turn to bare potentials bounded from below, but showing two minima. Presumably, this property also 
holds for the effective potential. However, since we employ a polynomial expansion of the potential we can only 
study the renormalization group flow around the electroweak minimum reliably. Future studies based on a different 
expansion will have to shed further light on the global renormalization group flow of the effective potential. 

As already discussed in Sect. |HCl even a very small positive value of Ag is sufficient to ensure stability. However, 
the point H = Q \s typically only meta-stable in this case. The phenomenological viability of this scenario requires 
that tunneling between the minima is sufficiently slow. To ensure this we enforce the conservative choice (cf. 
App.|^, 


\A{k) > -0.052, Vfc. 


(19) 


With this choice, the limit therefore essentially reduces to the usual longevity limit for the meta-stable region as 
obtained in the absence of any Ag. 

Our numerical study based on a more advanced approximation as described in App. 0 shows that already with 
moderately small values Ag « 0.25 the cutoff scale can be increased by approximately two orders of magnitude 
while retaining the full stability of the electroweak vacuum. This can be seen from the green region I in Fig. 
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A6(A) 


Figure 6: Different stability regions as a function of A and Xq. In region I (green) the potential is stable everywhere; in 
region II (blue) the UV-potential is stable, while the potential is only pseudo-stable for intermediate scales v < k < A-, 
in region III (red) the UV-potential violates the unique-minimum condition. Our simple model includes the effects of the 
electroweak gauge bosons only in the running of A4. In Appendix]^ we model the electroweak contributions to the full 
Higgs potential. The corresponding region III is then somewhat larger, cf. Fig. |12[ 


Increasing the possible cutoff scales by further orders of magnitude is difficult since there is a strongly infrared- 
attractive pseudo-fixed point at [Ael « 0 for those scales, see Sec. HE below. Next, there is a large pseudo-stable 
region (blue) where the UV-potential as well as the low-energy effective potential are stable, but our polynomially 
expanded potentials exhibit further minima at intermediate scales. As already explained, this is beyond the strict 
validity of our approximation and further studies are needed. Already relatively small values of Ag are sufficient 
to stabilize the UV-potential — although not the minimum at H = 0 — up to the Planck scale, since A 4 is 
negative, but its absolute value remains quite small. Finally, in the red region III the UV-potential is already 
meta-stable, as might be the effective potential. In this region the cutoff can easily take values beyond the Planck 
scale even for tiny values of Ag, because the tunneling rates are always small enough to guarantee the longevity 
of the electroweak vacuum. 

If we are well into region I, our RG flow indicates that the effective potential at fc « 0 is stable. Our approxima¬ 
tion should be reliable in this region. The fate of the pseudo-stable region II is not quite as clear. The appearance 
of the second minimum at intermediate values of k hints at the possibility that the effective potential at fc « 0 
might be meta-stable. However, to fully establish this feature requires an approximation that goes beyond the 
local polynomial expansion we have used. Finally, in region III it is likely that the effective fc « 0 potential will 
be meta-stable. 


E. Fixed-point structure 


To understand the behavior of higher-order couplings below the cutoff scale A, it is useful to analyze the beta 
functions in a little more detail. For couplings Ag and higher, they have the form 


/3a„ = (n - 4)A„ -b r]Xn + c , 


( 20 ) 


The first term reflects the canonical dimension of the operator. The second term is a generalized anomalous 
dimension term, i.e., r] is independent of A„. We also include contributions which are strictly speaking not an 
anomalous dimension, but are nevertheless linear in the coupling. For example, j5xg contains such a term scaling 
like A 4 Ag. Finally, c is a contribution from other couplings and hence independent of A„. For instance, top quark 
fluctuations contribute proportional to y". The beta function for Ag thus features terms at most linear in Ag. It 
is straightforward to find a pseudo-fixed point 


A 


6 * 


C 

2 -br/’ 


with 


2nyN^y^ - lOSAf 
IGtt^ 

QuyNctp + 90 A 4 


( 21 ) 


where the values for c and rj can be read off from Eq.([^. The pseudo-fixed point shows a remarkably small scale 
dependence, since the running of A 4 and y is also comparatively small. In fact, the ratio X^jy^ also approaches 
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Figure 7: Running Xeik) along RG trajectories with A6(feuv) = 0 (black solid), Xeikvv) = 0.1 (red dashed), A6(A:uv) = 0.5 
(blue dotted), A6(fcuv) = 2 (green dot-dashed) and the pseudo-fixed point (thin grey dashed). 


a fixed point [22] • We neglect a possible contribution from Ag. Because the anomalous dimension t] is small, 
the canonical dimension renders the fixed point strongly infrared attractive. In Fig. [^ we show that arbitrary 
starting values for Ag quickly converge toward the pseudo-fixed point. At low scales, Ag follows the rapidly growing 
top Yukawa coupling. As expected from the sign of the fermionic contribution in Eq. (HI), Ag becomes negative. 
This does not signal an instability, as higher-dimensional couplings also follow the value dictated by the Yukawa 
coupling and the fermion determinant contributes positively to the effective potential. The stability of the effective 
potential including top quark effects is thus guaranteed by our finite ultraviolet cutoff. 

In the perturbative approach, the explicit Ag dependence is not taken into account. However, some of the 
corresponding effects are encoded in higher loop diagrams [53]. This can be accounted for by setting Ag to its 
running pseudo-fixed point value [22j . In Fig. [^we also show the running fixed-point value from Eq.Q, i.e. the 
effective value of Ag included in loop effects in our definition of j3x^. To understand why the A„, n > 4 follow 
their respective pseudo-fixed point values on all scales in the perturbative approach, one should keep in mind that 
the cutoff is sent to infinity in this case. In an asymptotically free theory Ag would vanish in the far ultraviolet 
and correspond to a UV-repulsive direction. Therefore, to stay on an asymptotically free renormalization group 
trajectory would require a fine-tuned value of Ag in the infrared, such that one would lie exactly on the UV-critical 
surface. This corresponds to setting Ag essentially equal to the running fixed point value. Because the Standard 
Model is not asymptotically free, there is no reason to require that Ag should be set to its running pseudo-fixed 
point value in the ultraviolet. Toward the infrared, it will still approach the same value. In an effective theory Ag 
can be chosen independently of the pseudo-fixed point at A. A given ultraviolet completion will actually dictate 
the ultraviolet boundary condition. A similar statement applies to all higher-dimensional couplings. 


III. MODELS FOR HIGH-SCALE PHYSICS 


Starting with an effective held theory at a hnite scale A, the corresponding action can depend on many additional 
free parameters in the top-Higgs sector. An example is given by Ag as studied in Sec. [TT] These parameters can 
have a sizeable inhuence on the stability of the Higgs potential at high energy scales. In Sec. [HD we demonstrated 
that the interplay of just the Ag term with the Higgs quartic coupling shifts the limits on the measured Higgs 
mass from naive stability arguments by several per-cent. Let us now address the key question as to whether such 
an effect can be generated in a particle physics model, by integrating out heavy states. Depending on the choice 
of cutoff scale, we would of course expect that the microscopic model could also contain gravitational degrees of 
freedom. Here, we will focus on a much simpler toy model to demonstrate the presence of higher-order couplings 
at A. 


A. Heavy scalars 

As a first step, let us check if additional heavy scalars can provide a model for new physics above A, which 
features stable potentials with Ag(A) > 0 and A 4 (A) < 0. In the simplest scenario, the cutoff scale A corresponds 
to the mass scale of additional states, which are coupled to the Standard Model. For scales k > A they contribute 
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to the renormalization group flow of the SM couplings shown in Eq.Q. To compute the allowed Higgs masses 
in such a scenario, we do not need to consider a UV-complete theory for the heavy states beyond A. Instead, 
the model can come with an inherent cutoff, Absm 3> A. This corresponds to a hierarchy of effective theories, in 
which the Standard Model is superseded by a model containing heavy scalars, which again will be embedded in 
a more fundamental model close to the Planck scale. 


In our simple model we couple the additional scalar to the Standard Model through a Higgs portal 
to the effective potential of Eq. ^ , 


added 


AEeff(fc) = XHs{k)^Y + 


( 22 ) 


Due to the reflection symmetry S —)■ —S, the heavy scalar is stable and could constitute (a part of) the dark 
matter relic density. We assume that this Z 2 symmetry remains unbroken. The additional massive scalar field 
adds new loop terms to the beta functions for A 4 and Ag, contributing only for k > ms- To decouple the massive 
modes below ms we include threshold terms of the form 1/(1 + m'g/k'^)’^ with an appropriate power n. Including 
these threshold effects, the loop terms can be calculated in the ERG scheme, as shown in App. 0 If we allow for 
Ns mass-degenerate scalar fields with the same Higgs portal interaction, the one-loop beta functions of our toy 
model, Eq.Q, become 


^A4 


1 


-UyNcv'^ -I- 2 nyNcy^Xi -I- 9A4 


/^Ae — 2A6 -|- 


1 

IGtt^ 


2nyNcy^ + QuyN^y^XQ 


ro , 4 

— Ag -b C\gp -b 
108A|-b 90A4Ag 


^S>^HS 

4(1 -b m’g/k'^Y 

2(1 -b m'g/kX')^ 


(23) 


The Higgs portal contributions to the beta functions follow the usual pattern for bosonic fluctuations: receives 

a positive contribution, while for the sign is reversed, etc. Most relevant to our argument is the negative sign 
in /3;yg, because it implies a growing coupling toward the infrared, and thus a stabilizing effect. Looking at the 
effective potential we can generalize the stabilizing effect of heavy scalars beyond the leading terms in a polynomial 
expansion, namely 


d Vefi{k) ^ 1 

dlogfc 1 + XhsH^ 


(24) 


Through the overall positive sign heavy scalars indeed decrease the value of the effective potential toward the 
infrared. Eurthermore, the strength of the effect depends on the value of the Higgs field, and is largest at small 
field values of H. In the following, we will restrict ourselves to A 4 and Ag, but keep in mind that higher-dimensional 
couplings will be generated with alternating signs. While truncating this series at finite order could suggest either 
stability or instability, the contribution to the effective potential generated by heavy scalars remains stable, and 
features a minimum at vanishing field. 


Eor approximately constant Xnsik) the value of A 4 (A) will depend on — log(ABSM/A). Thus, the bare value 
A 4 (Absm) will determine A 4 (A). On the other hand, higher-dimensional couplings such as Ag reach values which 
are independent of Absm, as expected from their canonical dimensionality: As in the case of the pure Standard 
Model, the infrared value of Ag in the presence of a heavy scalar is determined by a pseudo-fixed point. In Fig. 
we demonstrate that the behavior of Ag is completely determined in terms of this strongly IR-attractive pseudo- 
fixed point, making the value of Ag(ABSM) irrelevant for weak-scale observables, as long as A ^ Absm- All three 
scenarios shown yield the correct values mh = 125 GeV and = 173 GeV. The value of Ag(A) only depends 

on the relevant and marginal couplings at that scale. As expected from universality arguments, Ag forgets the 
dynamics between Absm and A, different from the marginally relevant coupling A 4 , which does depend on the 
details of the dynamics at Absm- 

Comparing the full solution to the pseudo-hxed point trajectory in and beyond the Standard Model we indeed 
observe that setting Ag(A) = Aggjy; (and correspondingly for all higher-dimensional operators) is too restrictive. 
Instead, the value of Ag is determined by an interplay of two pseudo-fixed points: As long as k > A, Ag is determined 
by the pseudo-fixed point in the presence of new physics. Below the scale of new physics, Ag undergoes a transition 
to the pseudo-fixed point determined by the relevant and marginal Standard Model couplings. A caveat is that 
this analysis applies only to situations with small anomalous dimensions, i.e. in a perturbative regime. 


We can now determine the conditions under which A 4 (A) and Ag(A) assume values that according to Sec. HD 


yield Higgs masses below the conventional lower bound. First, we re-iterate our earlier observation that seeming 
instabilities in the potential in Fig. are artifacts of our truncation to A 4 and Ag. To fix A 4 and Ag we set 
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Figure 8: Running couplings in the presence of a Higgs portal coupling Xhs and Ns = 3. In the top row the boundary 
condition is A 4 (Absm) = 0, while in the second row it is replaced by a finite value. In the bottom row we illustrate a 
moderately big Xhs ~ 1 at the Planck scale, which is sufficient for an absolutely stable potential at all scales. In the right 
panels, we show the pseudo fixed point value for Ae in the presence of heavy scalars, X^ hs^ without heavy scalars, Ae. 


A 4 (Absm) such that the Higgs mass comes out correctly. We then solve the pseudo-fixed point equation for Ae, 
which determines its value at A. We obtain 


AeCA) 


-Uy^A) + 216Al(A) + NgX^Hsi^) 
4(167r2 -h V(A) -h45A4(A)) 


(25) 


A positive top Yukawa and a negative Higgs quartic coupling, which is the case of interest for the Higgs mass 
bounds, each reduce the pseudo-fixed point value. As the top Yukawa coupling grows toward the infrared, this 
negative contribution increases. Accordingly, a sufficiently large and positive value of Xhs is needed for a positive 
fixed-point value Ag, which is then depleted toward lower scales. In contrast, A 4 is driven to increasingly negative 
values by the fluctuations of the heavy scalar. This implies that while the generated value of Ag provides a 
potential that is bounded from below, it is not always large enough to yield a potential with a minimum at 
vanishing field values. The size of the quadratic Higgs term y? still decides whether the potential at the cutoff 
scale A is meta-stable, or features a global minimum at vanishing Higgs field. 

It is nevertheless possible to generate initial conditions corresponding to a stable bare potential with a minimum 
at a vanishing field value. As an example, we use A = GeV, where y{A) = 0.476 and A 4 (A) = —0.017 give 
a physical Higgs mass of 125 GeV. This value of A 4 (A) can be reached by adjusting A 4 (Absm)- To obtain a 
UV-stable potential with the measured Higgs mass value, we need 

NsXls ^ 24 , (26) 


for example corresponding to Xhs ^ 2 and 3 additional scalars. This large value of NsXhs needed to compensate 
for the factor 216 in front of the A 4 term in Eq.(25), which arises from the combinatorics of the respective scalar 
diagrams. Accordingly, a new physics scenario in which the new states are not combinatorically disfavored in 
comparison to the Higgs can accommodate larger values of Ag without needing large numbers of new states 
and/or large couplings. 
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In our example, the large value of the Higgs portal coupling implies that the heavy scalar sector will be driven 
toward a Landau pole not far above Absm- In such simple models one could therefore conclude that a UV-stable 
potential for the Higgs sector can be generated at the cost of a non-perturbative regime not far above Absm- 


B. Heavy fermions 

An obvious open question is the effect of heavy fermions. If their mass is generated through symmetry breaking 
at the electroweak scale, they will have a large Yukawa coupling. For example models featuring a heavy chiral 
fourth generation show a significantly reduced value of the possible cutoff scale A since the extra fermions would 
just add to the problematic effect of the top quark. Furthermore, such models are experimentally excluded through 
the Higgs coupling measurements at the LHC. 

We consider a model where the additional fermions are heavy, but their Yukawa coupling is small. In this setting 
we rely on an unspecified symmetry breaking mechanism at a high scale, affecting only the additional fermions. In 
our simple model we include such a mass term without discussing its possible origin. While this is straightforward 
in the present Z 2 model, in spite of the mass term explicitly breaking the Z 2 symmetry, an embedding of such a 
mechanism in the Standard Model with its SU (2 )l gauge invariance is less clear. The relevant Lagrangian terms 
for Njj heavy fermions then take the form 


^ VnHTjTj + mlijr] . (27) 

Alternatively, we could consider a model where the additional fermions are singlets under the Standard Model 
gauge groups. In that case they can have a mass even in the unbroken phase. Their Higgs couplings would 
correspond to a dimension-5 operator of the form H^fjr]. 


Let us again analyze the induced potential at A 
of Eq.(25) now gets the additional contribution 


in terms of the pseudo-fixed point for Ag, which instead 


A6(A) 


-12^/6(A) + 216Al(A) + NsX^sW " 12IV^y^(A) 
4 (I67r2 -h V(A) -h 45A4(A) -h 9N^y^(A)) 


(28) 


The fermion 77 induces a negative contribution to Ag. In general, its contribution to A„ will be positive for n/2 
even, and negative for n/2 odd. This means that heavy fermions will generically make it hard to reach sizeable 
positive values of Ag, but they might be of interest for cases where, e.g., Ag > 0 stabilizes the potential. 


Alternatively, heavy fermions can be important in a scale invariant theory in the far ultraviolet, H(Absm) = 0. 
Here, fermions help to control the induced value of A 4 (A). In a simple analysis the masses of heavy bosons and 
fermions split the new physics scale A into two thresholds ms = A 5 and = A,j, where we assume ms < mr^. 
The induced value of A4 for Ag < A^ now reads 

A 4 (As) = A 4 (Absm) + log + ^y^(A^) log , (29) 

where we assume a slow running in the Yukawa and the portal couplings to the heavy scalars. A hierarchy 
> ms implies that the value of the induced coupling Ag will be independent of y,, and instead be determined 
by the pseudo-fixed point at Ag. In this framework we can construct models where A4 is small and negative and 
Ag = 0{1), without the need for large fermion and boson numbers or large couplings. 


C. Non-perturbative new physics 

In the last two sections we have shown that simple models with additional heavy fermions and bosons allow 
us to generate initial conditions for the renormalization group flow which yield low Higgs masses based on a 
UV-stable potential. The main challenge is to avoid the non-perturbative regime beyond Absm while generating 
coupling values jAgl > IA 4 I. In general, the effect of new physics on the beta functions of the running couplings in 
the Higgs sector follows two different patterns: 

1. terms independent of the Higgs self-couplings induce these couplings, even if they vanish at some scale Absm; 

2. terms which change the scaling of the Higgs self-couplings by inducing an anomalous dimension t/np- 
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Figure 9: We show A 4 (black) and Ae (purple), obtained using the model beta functions from Eq.( |30[ ), which include 
new-physics effects beyond A = GeV (dotted vertical line), induce these couplings and alter their canonical scaling 
dimensionality. 


This means we can write the beta functions for A 4 and As in the presence of new heavy states as 


87r2 
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—UyNcu'^ + 2 nyNc.y ^+ 9A4 —^4 
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/^Ae — 2 A 6 + ^^^2 — IO8A4 + 9OA4A6 + CgNp] + %NpA 6 , 


(30) 


where C 4 / 6 NP and ? 74 /gNp encode the effects of the microscopic degrees of freedom beyond the Standard Model. 
Again, it is useful to analyze the running of Ag in terms of a pseudo-fixed point of the kind shown in Eq.(21), 
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2 - 1-77 
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2nyNcy^ - 108A| -I- cgNP 
IGtt^ 
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+ 776 NP • 


(31) 


For example, we can start with vanishing values for all couplings at Absm- We can induce a negative value of A 4 
with increasing magnitude and a positive value of Ag, if we choose C 4 np >0 and cgNP < 0 and |cgNp| > ‘^UyNcy^. 
The pseudo-fixed point moves toward smaller values as A 4 grows in magnitude toward the infrared. As long as 
776 NP is small the pseudo-fixed point stays strongly infrared attractive. 

Let us now broaden our approach and allow for new physics in a non-perturbative regime, such that large 
anomalous dimensions of either sign can exist. An example could be asymptotic safety [39]. We choose a scenario 
where 77gNP — —2, i.e. , the scaling dimension of the higher-dimensional coupling changes from irrelevant to 
marginal, or even relevant. In such a case, the pseudo-fixed point will not be strongly infrared attractive anymore. 
Accordingly, new physics can induce a sizeable value of Ag, that is not depleted between Absm and A. From Fq.(31 1 


it is clear that cancelling the canonical dimensionality by a large anomalous dimension can lead to significantly 
enhanced values for Ag. 

As an example, we set C 4 np = 9/4, cgNP = —27/2, 774 NP = 0, and 77gNP = —2 and assume that the effect of new 
physics appears beyond A = 10^“^ GeV. We then obtain A 4 (A) = —0.026 and Ag(A) = 0.51, as shown in Fig. 
This does provide a UV-stable potential with Higgs masses signihcantly below the conventional lower bound. 


IV. CONCLUSIONS 

With the discovery of the Higgs boson the Standard Model describing the interactions between fundamental 
states in terms of a gauge theory is complete. The crucial question becomes whether the Standard Model itself 
gives us a hint where new physics must appear |3D]. One possible clue is the behavior of the Higgs potential at high 
energy scales and large field values. A perturbative analysis of the Higgs potential in the Standard Model indicates 
that it loses its stability around 10^° GeV, well below the Planck scale. We have investigated the robustness of 
this claim in the presence of higher-dimensional operators in the effective Higgs potential, as one would expect in 
theories with a physical UV-cutoff. Because of their canonical dimension the effects of these higher-dimensional 
operators will vanish quickly toward the infrared, but nevertheless can play a crucial role in the vicinity of the 
cutoff. We identify several scenarios: 
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• The UV-potential including higher-dimensional couplings with A 4 < 0 is stable, as is the effective potential. 
The electroweak minimum remains the global minimum throughout the RG evolution. In this case, already 
the operator allows us to extend the UV-cutoff by about two orders of magnitude. Further higher¬ 
dimensional operators may increase the cutoff scale further. 

• The UV-potential and the effective potential are stable with their respective global electroweak minima. 
However, during the RG evolution this minimum seems only meta-stable for an intermediate range of scales. 
For conclusive statements about this pseudo-stable scenario, a more complete treatment of the RG evolution 
of the effective potential is necessary. We expect that part of this parameter space belongs to the class of 
fully stable scenarios whereas other parts may feature a real meta-stability. Taking the scenario literally, 
the barrier separating the two minima of the fc-dependent potential is typically of the order k. 

• The UV-potential is meta-stable, featuring a global minimum at some high field value. This can be achieved 
for an arbitrarily small initial Ag, and therefore no instability occurs even for arbitrarily large negative A 4 . 
Implementing the constraint that the electroweak minimum has to be sufficiently long-lived, this essentially 
reduces to the standard discussion of a viably long-lived meta-stable region. 

Higher-dimensional operators affect the constraints on the Higgs mass between ^ 1 GeV for the stable potential 
to ^ 4 GeV for a pseudo-stable setup. In the presence of an ultraviolet cutoff, the size of higher-dimensional 
couplings is an inherent ambiguity of the Standard Model. To include higher-dimensional operators in our setup 
we have used non-perturbative functional renormalization group techniques in the presence of a finite UV-cutoff. 
Neglecting higher-dimensional operators, our approach and the usual perturbative approach agree in the domain 
where perturbation theory with a massless regularization scheme is applicable, i.e. , for energies and field values 
much smaller than the UV-cutoff A. With regard to higher-dimensional operators, perturbation theory assumes 
very specific and small values, essentially those generated by the Standard Model degrees of freedom and encoded 
in higher-order effects. Our non-perturbative treatment allows for more general values of the higher-dimensional 
couplings. 

The difference between the values of the couplings in perturbation theory and our more general values is 
in principle measurable. These effects become sizeable at energies where the new dimensionless couplings are of 
order one. Restricting the higher-dimensional couplings to be of order one at the cutoff scale (where they typically 
have their largest values), limits the size of the achievable shifts, but even in the absence of higher-dimensional 
operators, stability up to the Planck scale is only disfavored by ~ 3cr [3]. Small shifts as observed in our analysis 
can significantly reduce this tension. 

Physics beyond the Standard Model predicts the size of the Higgs self-interactions at the cutoff. We can ask 
what type of new physics generates higher-dimensional operators of suitable size, such that Higgs masses below the 
conventional lower bound can be obtained. In general, UV-completions featuring large anomalous dimensions can 
easily generate sizeable higher-dimensional operators, stabilizing the Higgs potential up to the Planck scale. Their 
effects can be linked directly to the more general analysis in terms of higher-dimensional operators contributing 
to the Higgs potential. 

In weakly coupled models the answer is less obvious. Our approach to new physics affecting vacuum stability 
is different from the usual strategy of modifying the running of the Higgs potential through new particles below 
10^*^ GeV, as discussed in the Appendix. In the main body of our paper we focus on very heavy new particles, 
coupled to the Standard Model through renormalizable operators. They modify the Higgs potential through 
induced higher-dimensional Higgs self-couplings. Specifically, we have investigated a simple extension of the 
Standard Model by heavy scalars and fermions. While it is possible to generate sizeable higher-dimensional 
operators, the required parameter choices in our simple model (couplings and number of fields) are on the border 
of being non-perturbative. It remains an intriguing question, what kind of new physics can generate higher¬ 
dimensional operators of a suitable size, while remaining perturbative on all scales. 
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Appendix A: Light dark matter scalars 


Different scenarios predicting new particles which modify the relation TOij(A) from the conventional lower 
bound can be distinguished by the ratio of their mass to the cutoff A. In the main body of this paper we focus on 
heavy states with mass m> A. Their fluctuations induce an effective potential 14fj(A) which generically contains 
non-vanishing higher-dimensional terms. 

A different scenario arises using light states with mass m > v. Their primary effect is to modify the running 
of (some of) the Standard Model couplings between the electroweak scale and the cutoff scale. Most importantly, 
they can alter the RG flow of the quartic Higgs self-coupling, such that it remains positive all the way to the 
cutoff A. One example is a gauge singlet, dark matter scalar. It is coupled to the Standard Model through the 
Higgs portal defined in Eq.(22), but now with a mass ms in the GeV-TeV region [TT1[24]. With a Z 2 

symmetry protecting it from decaying and in the absence of a second VEV which would lead to mixing with the 
Higgs scalar, the only additional parameters of this model are ms and Xhs- This makes it the arguably simplest 
model, which can accommodate the complete observed dark matter relic density |25j . Erom a field theory point of 
view, the portal coupling cannot be large and negative to avoid an unstable combined scalar potential. To yield 
the desired effect on the Higgs mass bound it should be positive. 


The difference to the additional heavy scalar discussed in Sec |HI A] is that Xhs is now important at scales k < A. 
A sufficiently large Xhs can altogether prevent A 4 from becoming negative. In Fig. 10 we update results from 
Ref. m by including the effect of the strong gauge coupling and approximating the effect of the weak couplings 
through a fiducial contribution as defined in Eq.Q. A moderately large Xhs ~ 0.78 at the Planck scale is indeed 
sufficient for an absolutely stable potential for the Higgs at all scales. In our updated analysis we also guarantee 
the correct weak-scale masses mh = 125 GeV and = 173 GeV. The infrared value of Xhs in Fig. 


10 


IS 

compatible with a dark matter mass of ms = 340 GeV, if we require the new scalar to constitute the complete 
dark matter relic density [25j . 

From this example we tentatively conclude that the positive contribution of Xhs to is sufficient to reconcile 
the measured value of the Higgs mass with a stable Higgs potential on all scales up to the Planck scale. As a 
crucial phenomenological distinction to our investigation in the main part of this paper, new physics states appear 
at low scales ms ^ v, instead of at mass scales ms ~ A. Thus the new light states can become experimentally 
accessible in dark matter searches of at colliders in the near future. The particular region of parameter space used 
in Fig. 10 will be probed by planned upgrades of XENON [2S] and LUX |27) . 


Appendix B: Tunnelling 


In this appendix we briefly discuss estimates for the tunneling rates in the presence of a ^ Agi/J^-term [5T]. The 
tunnelling rate can be calculated by determining the Euclidean action of bubbles S'Elv^bounce]) where (/^bounce is 
the bounce solution of the Euclidean equations of motion for a real scalar field 
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ounce 
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Figure 10: We show the running of the couplings employing Eq.Q, amended by an additional Higgs portal coupling Xhs- 
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with the boundary conditions <;3bounce(oo) = <i3bounce = and V^bounce(O) = 0- The value is the location of 
the meta-stable minimum. The tunneling rate then becomes 


r 

V 


ounce] 


Let us now turn to the potential given by Eq. namely 


V((p) 


ml 






(B2) 


(B3) 


For the formation of a Bubble of size ^ R one should evaluate the scale dependent potential at the scale k ^ 
1/R pSl m]. In the situation we are interested in, the field values where the potential starts to be lower than 
the electroweak vacuum near zero vev are quite large. Accordingly we are interested in very small bubbles and 
k ^ 10^° GeV. Neglecting higher-order terms, the potential would be unstable with ml > 0 and A 4 < 0. To 
stabilize the potential we can add a term with, Ag > 0 . 


We are mainly interested in an upper limit on the tunneling rate. It therefore makes sense to approximate the 
potential Eq. (B3) by a potential which is always lower than the potential we actually have. A simple choice is to 
simply take the term 


Vapprox(</?)= -A4<P < E(v3). (B4) 

For such a potential the tunneling rate has been computed in Ref. [23]. In this case the bounce solution is. 


V^bounce 
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Because the approximate potential is scale invariant we actually have infinitely many bounces corresponding to 
different values of i?, and in principle we have to take all of them into account. 

To determine the pre-factor of the exponential one has to perform a one-loop calculation around the background 
of the bounce solution [3D|. Doing this |23| one finds for the tunnelling probability in a space time volume V 4 


V4 / Stt^ a 

RA®''P'v"3|A4(l/R)| j ■ 


(B 6 ) 


The scale dependence of A 4 breaks the scale invariance and differently sized bubbles have different actions. To de¬ 
termine the loop-correction one has to calculate the inverse of the determinant over the spectrum of the fluctuation 
modes around the bounce background. The zero-modes have to be treated with some extra care. They correspond 
to integrations over the corresponding collective coordinates. The volume factor arises from the translation zero- 
modes and the corresponding integration over the space-time volume. The scale invariance of the approximate 
potential leads to an additional zero-mode in the spectrum. The integral over dR/R is the integration over the 
collective coordinate corresponding to this zero-mode m- Taking all possible bubbles into account we finally 
arrive at 
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dR V 4 f Stt^ a 
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(B7) 


For constant A 4 (l/i?) this integral is divergent for small values of R. Such small values of R correspond to a large 
maximal field values of the bounce solution, as can be seen from Eq.(B5|. Indeed we have. 


V^bo 




(B8) 


This tells us how the integral is regularized in presence of Ag > 0 which stabilizes the potential. For field values 
beyond the point where the potential has its minimum and starts to rise again, the approximation of the potential 
by just the —A 4(/?^/4 term is not sensible anymore and we should cut off the integral at this point at the latest. 

In practice, the integral is typically dominated by a reasonably sized region around the most negative values of 
A 4 because of the strong exponential suppression of the integrand for small IA 4 I in Eq.(B7|. Even a small upturn 
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in A 4 (l/i?) typically quickly overcomes any growth in the 1/ii^ factor. Indeed a better approximation is probably 
to use an effective quartic coupling 


A 4 ® - ^4 + (B9) 

Here, one can again see the stabilizing effect of a positive Xq. With increasing field value, i.e. decreasing R of the 
bounce, A 4 ® becomes less negative effectively cutting off the integral. More intuitively, an increasing Ae lowers the 
depth of the minimum, and thus lowers the tunneling probability. 

Let us now determine the value of the meta-stability bound on A 4 . For the space-time volume we can take the 
Hubble volume and time. 


^4 ~ tL, (BIO) 

where Hq is today’s value of the Hubble constant. Longevity of the vacuum requires that the tunneling probability 
in the volume V 4 ~ ^/Hq is smaller than 1. For R between R = 1/(10^° GeV) and R — 1/(10^® GeV) this 
corresponds to a bounce action Se ^ 400 — 500. Therefore, as long as 

|A4|< 0.052 (Bll) 

we are on the safe side. In the presence of higher-dimensional operators, this inequality will be relaxed, and thus 
represents a very conservative upper limit. 


A similar argument can be made when the higher dimensional operator is actually destabilizing, i.e. , Ae < 0. 
In this a suitable approximate potential is 


Amin \ 

Vapprox(‘/5) = with Ar" = ^4 + ^ . (B12) 

where (/Smax is the value where the potential is once again stabilized by even higher dimensional terms. Using the 
same arguments as above the limit becomes 
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< 0.052. 


(B13) 


We conclude, in accordance with Ref. [51] that the presence of higher-order operators could have a significant 
impact on the lifetime of the electroweak vacuum in a meta-stable situation, for example if Ag is negative. Estimates 
relying solely on A 4 strictly only apply to a restricted set of UV-completions. Large higher dimensional operators, 
in particular ones that destabilize the potential can have significant effects. 


Appendix C: Functional renormalization group 

The /3-functions for the running couplings used in this work follow naturally from the functional renormalization 
group (RG). Formulated in terms of a flow equation for a momentum-scale-dependent effective action Ffc, the 
functional RG provides for efficient means to follow the flow of higher-dimensional operators, to account for 
threshold effects and to even approach regimes of strong coupling. Of course, the universal one-loop /3-functions 
dominantly used in the main text follow as a simple by-product in the weak-coupling limit. 

In the following, we briefly recall the essentials of the functional RG, concentrating on its application to the 
present model. As this method has already been employed in the present context to the pure Z 2 Higgs-Yukawa 
model m, we focus on the generalizations for the SU{3) gauged version. Finally, we detail how the universal 
one-loop /3-functions arise by neglecting higher-loop terms as well as RG improvement and by considering the 
‘deep Euclidean region’. 

The functional RG is a manifestation of the Wilsonian idea of solving a quantum field theory by integrating 
out fluctuations in the path integral, momentum shell by momentum shell. In its modern version, it can be 
conveniently formulated in terms of a scale-dependent action functional Ffc. This is constructed such that it 
interpolates between the UV-theory parameterized at the cutoff scale A in terms of the Ffc_>.A —>■ Sa and the full 
effective action in the IR, Ffc_>o —>■ Feff. The latter contains, for instance, the full effective potential as a local term 
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without derivatives of the fields, Feff = / + 0((d</>y)j. The change of the effective action as a function 

of the RG scale k (RG flow) is governed by the Wetterich equation [18] 


dtTk = -STr 

t fc 2 


dtRk 
+ Rk 


t = In 


A 


(Cl) 


where denotes the second derivative of with respect to the fluctuating fields {H, ijj, ip, A^, •••); the 
super-trace includes a minus sign for fermionic trace parts. It also involves a summation over the eigenvalues of 
the regularized propagator, -\- Rk)~^, which in our case reduces to a momentum integration. A new technical 
ingredient is given by the regulator Rk which can be thought of as a momentum-dependent and fc-dependent 
mass term. Specifying Rk corresponds to specifying the details of the momentum shells regularization. Different 
choices of Rk correspond to different regularization schemes. 

Even though the Wetterich equation has a simple one-loop structure, it is an exact equation as there exists an 
exact equivalence to the full functional integral. The key ingredient for this is that the propagator in the loop 
^ -I- Rk) denotes the full propagator at the scale k. From the full solution of Eq.(Cl), for instance, all 

correlation functions can be computed. 


While a perturbative expansion of the Wetterich equation reproduces perturbation theory to any order, the flow 
equation can also be used to extract non-perturbative information m- Following the Wilsonian idea, the scale 
dependent action Fk contains not just relevant or marginal couplings, but encodes effects of high-scale quantum 
fluctuations in the presence of higher-order operators. For practical calculations, this infinite series of operators 
in the effective action is truncated. In the present context, a useful ordering scheme defining our approximation of 
the exact flow is given by an expansion of the effective action in terms of derivatives of the field. To next-to-leading 
order this gives 


Fk = 




rr 

V+^{d^H)^ + ^Z^^iP,tlpiP, 




(G2) 


where the potential V, the (bare) Yukawa couplings yj as well as all wave function renormalizations Zh,^,g 
are fc-dependent. The gauge-part of the action is also supplemented by a gauge-fixing and Fadeev-Popov ghost 
contribution. 

Inserting this ansatz into the Wetterich equation and expanding both sides in terms of this basis of operators 
leads to the /3-functions, i.e. the flow equations for the Yukawa couplings and the wave function renormalizations. 
For the flow of the effective potential, we obtain a /3-functional defining a ,5-function for every field value H. For 
clarity, let us write down this flow of the effective potential explicitly: 
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Here, we have introduced the anomalous dimensions of the fluctuation fields, defined by the wave function renor¬ 
malization flow. 
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In Eq.(C3), we have also made an explicit choice for the regulator function in terms of the linear regulator |42] . 
Other choices lead to a qualitatively similar behavior on the right hand side. This expression clearly contains 
threshold effects, where for example the scalar contribution will be dampened once the scalar mass, related to 
V”, become s of order k. 

From Eq.(G3l, the /3-functions for the scalar self-interactions X 2 n can be obtained straightforwardly. To do this, 
we expand the potential in powers of the field. 


VMk) = 2 
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and introduce the renormalized couplings. 
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Here, the rescaling with the wave function renormalizations compensate for the field rescalings which are necessary 
to bring Eq.(C2) to a canonical form, e.g., H —^ H/Z^ , etc. Inserting the expansion of Eq.(C5l into Eq.(C3), 


and expanding both sides to fourth order in the field, we can read off the flow of A 4 , 


/3a 4 = 2?7_f/A4 + 
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We observe that the scalar loop terms ~ A| and ^ Xq are suppressed if the mass parameter is larger than the 
RG scale k. This signals a typical threshold effect characterizing the decoupling of massive modes. In the present 
application, however, is chosen such that it remains small for all values of k and eventually drops below zero 
near the weak scale /cew- Then all modes become massive and decouple, which can be made explicit in the flow 
by expanding about H = v. 

For the high-energy applications discussed in the main text, we take the limit of the deep Euclidean region, 
—>■ 0, as is standard in perturbative computations. Using the one-loop result for the anomalous dimension 
r]H (also obtained from the flow equation by studying the flow of the scalar kinetic term), t]h = -^c 2 /|/( 47 r)^, 

we obtain the I3\^ function as given in Eq.Q, except for the last term ~ gp modelling the electroweak gauge sector 
which is described in the main text. In fact, the ansatz of Eq.(C2) also yields higher-order terms like ^ X^rjH, 


corresponding to contributions from higher loops. The fact that these terms arise from the one-loop formula in 
Eq.( |Cl[ ) is due to the ‘RG-improvement’ of the propagator in the loop. 

For the perturbative estimates discussed in the main text, these higher-loop contributions are neglected, cf. 
Eq.Q, whereas they are fully included in the ERG results of F ig. [ 6 | further outlined below. 

Similarly, of Eq.Q is obtained from an expansion of Eq.( |G3[ ) to order H^, taking the deep-Euclidean limit 
and ignoring RG-improvement terms. Also a contribution ~ Ag is ignored in Eq.Q, whereas it is included in the 
ERG flow studies, see below. The derivation of Py of Eq.Q with yj —>■ y for Uy = 1 follows the same pattern. 
In the general case, the flow of the Yukawa couplings for different fermions can be derived in this way; then, 
the different fermions also acquire separate anomalous dimensions. The standard one-loop j3g^ function for the 
gauge sector can most easily be derived from the flow equation, using the background-gauge formalism. Here, the 
running of Qs is directly related to the anomalous dimension, j5g^ = rjGgs/‘^ see e.g. [H]. 

The manner by which we extracted the one-loop /3-functions from the Wetterich equation follows a general 
pattern: first expand both sides in terms of the operators under consideration; this yields the generally non- 
perturbative /3-function of the corresponding coupling within a given ansatz. For the one-loop order, all RG- 
improvement terms corresponding to higher-loop corrections can simply be dropped. The universal one-loop 
coefficients arise in the limit of the deep Euclidean region which is the expected necessary prerequisite for one- 
loop universality. While our explicit examples discussed above have used the linear regulator as a special choice, it 
can be proved that the one-loop coefficients for the renormalizable operators are manifestly scheme independent. 
However, the /3-function coefficients for Ae do depend on the choice of the regulator. This scheme dependence is 
physical in the sense that it parameterizes how the Standard Model as an effective theory is embedded into an 
underlying theory. 

For the non-perturbative functional RG results in the main text we integrate the /3-functions for the full ansatz 
Eq.( |G2[ ), including the mass parameter y^{k), all threshold dependences and all RG-improvement terms due to 
the back-feeding of anomalous dimensions. For the flow of the effective potential we use a polynomial expansion 
to order iJ®. Higher orders can easily be included, but do not lead to any significant change of results. For the 
pure Z 2 model, the local convergence has been checked to order ~ X4oH*° M- 

The flow equations for the pure Z 2 model are given in Refs. mis] for a general regulator function and the 
linear regulator. To match them to the present matter content, all fermion loop contributions have to be multiplied 
by a factor -^c, corresponding to the number of fermions interacting with the scalar sector. In the following, 

we confine ourselves to listing only the new contributions induced by the SU{Nc.) gauge sector. The electroweak 
sector will be modelled in terms of a fiducial coupling, as in Eq.( |7t. 

Whereas the form of the flow of the effective potential in Eq.(G3) or Eq.(12) in Ref. [T3] remains unaffected, 
the flow of the Yukawa coupling receives direct contributions from the gauge sector. Also for numerical purposes, 
it is convenient to write this /3-function as 


dy: 


d log/c 


= 2xEq.(13) of [H]- 
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The factor of 2 in front of the first term is due to different conventions for the Yukawa coupling. The dictionary 
for the conventions used in this work and in Ref. [13] reads: y^=2h‘^, /ZH=2k‘^p, V{H)=k'^u{p), where the 

LHS corresponds to this work and the RHS to M- The dimensionless quantity k is related to the expansion point 
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of the potential in field space; in the present work, we have k = 0 aX high-energy scales and k = v‘^/{2k'^) in the 
broken regime near the weak scale. The threshold functions I (and m below) describe the (regulator-dependent) 
decoupling of massive modes; for any admissible regulator, they approach finite constants for vanishing first 
arguments and decrease to zero for large first arguments. The result holds for arbitrary Lorenz-gauge parameter 
All numerical studies have been performed in the Landau gauge 5 = 0. 

Correspondingly, the fermion anomalous dimension reads 


Vi,, = Eq.(16) of [U 
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while the corresponding equation for the scalar anomalous dimension r]H remains form invariant to Eq.(15) of 

Ref. [H]. 

In practice, the resulting system of coupled differential equations is solved numerically. For stability reasons 
in the presence of higher-dimensional operator, the flow is solved from the UV down to the weak scale. The 
boundary conditions in the UV are chosen implicitly such that the infrared observables are matched as described 
below. 


Appendix D: Mapping between flow trajectories 


To study the RG evolution of the Standard Model, its physical masses and couplings have to be related to the 
fundamental parameters appearing in the /3-functions. Within the perturbative setting this procedure is well- 
established up to NNLO [5]. For our perturbative results we follow these prescriptions at the corresponding loop 
order, setting the parameters at the universal RG scale k = such that they reproduce a running top mass 

rrit = 164 GeV, a running Higgs mass rrih = 125 GeV, and the strong coupling as{mz) = 0.1184. We list the 
explicit numbers employed in our calculations for our model with Uy = 1 in Tab. |T] 


The functional RG approach allows us to access the regime of spontaneously broken symmetry, where a non¬ 
vanishing minimum iLmin.fc 0 oi the effective potential appears in the RG flow. This way all particle masses are 
generated dynamically. In contrast to the perturbative setting, the running minimum enters the RG evolution of 
the Higgs quartic coupling and the Yukawa coupling explicitly: The generated running masses mt(k) cx yH,ain,k 
and mh(k) oc y^ffmin.fc induce a threshold behavior that leads to a freeze-out of model runnings in the infrared. 


as shown in Fig. 11 The ERG flow for the RG relevant direction parameterized by the mass parameter has to 


be fine-tuned at the UV-cutoff scale, such that the minimum of the effective potential approaches the SM value 
of the vev in the infrared iLmin,fc->.o = v. This fine-tuning is a manifestation of the hierarchy problem, which is 
dealt with here by an explicit choice of UV-parameters. 

The threshold behavior in the ERG is non-universal, which means it depends on the choice of regulator function. 
This dependence is natural as it parametrizes the details of decoupling of massive modes from the flow. In the 
symmetry-broken regime, the threshold behavior leads to a significant deviation of the running couplings from 
their values in the perturbative setting, which is strictly tied to the symmetric deep Euclidean region; the latter 
is, of course, inappropriate for the threshold region. To compare to the perturbative results and extract universal 
physical results, we have to establish a matching procedure. We choose to match the flows at the scale 1 TeV, 
which is sufficiently deep inside the Euclidean region but not yet affected by possible higher-dimensional operators, 
cf. TabU _ 

Our procedure equates the running couplings in the perturbative MS scheme and the ERG scheme at the 
matching scale. Approaching that scale from the ultraviolet, typical ERG trajectories are still in the symmetric 
regime, an electroweak minimum of the potential has not been generated yet, and the threshold effects in the 
ERG /3-functions are subleading. In Fig. El we see that this regime features flows of the Higgs self-coupling and 
the top Yukawa coupling closely resembling the perturbative RG flows. As the values for the couplings agree with 
those of the MS scheme at the matching scale, our matching scheme guarantees that UV-initial conditions for the 
ERG flow are mapped to physical values in the infrared with high accuracy. Qualitatively, this mapping works 
already at the level of our ERG trajectories without a matching procedure, since threshold effects automatically 



Ar y Qs 

k = 

0.129 0.940 1.167 

k^l TeV 

0.101 0.867 1.060 


Table I: Values of the model parameters at the top pole mass and the RG matching scale k = 1 TeV for the matching 
between the perturbative RG and the functional RG trajectories. 
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Figure 11: Comparison and matching of FRG flows with the perturbative setting. The left panel distinguishes the symmetric 
regime and the regime of spontaneously broken symmetry. The center and right panels show the FRG flows (red dotted) 
as well as the perturbative flows (black curve) of the Higgs coupling and the top Yukawa coupling. The running couplings 
are matched at 1 TeV. 


provide for a freeze-out of the RG flow, such that in principle physical values can be extracted at fc = 0. In a 
more elaborate FRG setup, using a more sophisticated truncation of the operator space, our mapping using MS 
parameters will become obsolete and quantitatively precise values for physical observables can be read off from 
the FRG flow trajectories at fc = 0 directly. 

Our choice of matching scale is also applicable when the symmetry-breaking scale lies at around 1 TeV, as 
threshold effects remain small close to that scale. Matching the flows at 1 TeV therefore provides us with a 
well-defined procedure to compare the two RG schemes. As an added benefit we can thoroughly study the FRG 
evolution of our toy model consistently including higher-dimensional operators and threshold effects in a Standard 
Model context. 


Appendix E: Modelling electroweak contributions to the full effective potential 


In the present work, the electroweak contributions are modelled by means of a hducial gauge coupling con¬ 
tributing to the perturbative /3-functions Eq. 0 . In particular, we neglected a similar contribution to the flow 
of Ae in this simple model. With the full /3-functional for the full effective potential Eq. (C3) at hand, we can 


in fact model the electroweak gauge contributions to the full potential including the flow of the mass-parameter 
H^{k) as well as all higher couplings A 6 , 8 ,...(fc). To lowest order, we add the contribution of a hducial gauge loop 
to the /3-functional [35], 


Eq. (G3I 


Eq. (G3I 


cv 


327r2 


1-b 


4k^ 


(El) 


with the hducial gauge coupling gp and a numerical constant cy. As in the main text, we choose these hducial 
parameters in such a way that the running of the quartic coupling A 4 is quantitatively similar to the running of the 
Standard Model Higgs quartic coupling, cf. Eq. Q and Fig.[^ Comparing to Eq. Q, we identify cy = 16ca = 9. 
This choice leads to the same how equation for the quartic coupling as in the main text. In addition, also the how 
equations of all other couplings in the scalar sector {p?{k), A 6 , 8 ,...(fc)) receive contributions proportional to an 
appropriate power of the hducial gauge coupling gp. The flow of all scalar couplings is only slightly modified such 
that the solution for the simple model lies almost on top of the results with the advanced hducial modihcation 
Eq. (El). The only relevant quantitative change is the flow of the mass parameter /i^(fc). The additional bosonic 
degrees of freedom reduce the ratio g?{k)/k'^ on a wide range of sc ales . 


The main modihcation introduced by the advanced model Eq. (El) affects the validity regime of the unique- 
minimum condition Eq. (14). This leads to a shift of the border between region II and III of Fig. Smaller values 
for /r^(A) have to be compensated by a smaller UV cutoff A for hxed A 4 (A) and A 6 (A), resulting in a depletion of 
the transition curve. This effect is less pronounced at the transition between region I and II because the relative 
reduction in g?/k"^ is smaller. 
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Figure 12: Different stability regions as a fnnction of A and Ae analogous to Fig. [^for the simple model. In region I (green) 
the potential is stable everywhere; in region II (blue) the UV-potential is stable, while the potential is only pseudo-stable 
for intermediate scales u < fc < A; in region III (red) the UV-potential violates the unique-minimum condition. The 
inclusion of electroweak contribntions to the full potential leads only to minor modifications of the border between region I 
and II. Region III becomes somewhat larger as the smaller mass parameter facilitaes a violation of the unique-minimum 
condition already for smaller cutoffs. 
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